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УДК 517.9--528.21--550.831 
А.В. ЧОРНИЙ, д-р фіз.-мат. наук, ін-т геофізики НАН України 


ПРО НОВУ ЗАДАЧУ ДЛЯ РІВНЯННЯ ЛАПЛАСА 
Розглядається задача відновлення потенціалу притягання за значеннями модуля градієнта потенціалу сили 
тяжіння. Її розв'язок подається у вигляді послідовності розв'язання зовнішніх граничних задач Неймана для 
рівняння Лапласа. Доведена збіжність послідовності до розв'язку, що розшукується, за умови, що він не дуже 
відхиляється від заданого. 


Для розв'язання важливих теоретичних і прикладних задач різних наук про Землю, що виникають при 
вивченні її фігури, внутрішньої будови, глибинних геофізичних процесів тощо, а також задач, пов'язаних із 
рухом штучних апаратів у навколоземному просторі, потрібно знати розподіл значень потенціалу сили 
тяжіння, значень градієнта цього потенціалу чи значень модуля цього градієнта. Такий розподіл можна 
одержати на основі даних спостережень якої-небудь характеристики поля 1 фундаментальних властивостей 
функцій, що описують гравітаційну взаємодію. 

У цей час на поверхні Землі вимірюються в масовій кількості значення сили тяжіння, які являють 
собою значення модуля градієнта гравітаційного потенціалу. Однак до цієї пори ще не створена теорія від- 
новлення потенціалу за значенням його градієнта, У статті розробляється один із можливих способів роз- 
в'язання цієї задачі, який реалізується за умови, що відновлюваний потенціал не дуже відхиляється від 
заданого. 

Вступ до теорії задачі відновлення потенціалу за модулем його градієнта розпочнемо з описання 


солютно тверде тіло, за формою близьке до тіла обертання, що рухається рівномірно вздовж своєї орбіти і 
обертається навколо своєї осі з постійною кутовою швидкістю (без прецесії і нутації). Позначимо через у 
обмежену область точок тривимірного евклідового стору, яка зайнята масами Землі (за винятком мас її не- 
рухомої газової оболонки), через у | - необмежене доповнення цієї області, вільне від будь-яких гравітую- 
чих об'єктів, а через ду - границю між множинами у |і у , що ототожнюється з фізичною поверхнею 
Землі. Міри кожної з підмножин певної з областей у чи у відповідно до наданої моделі залишаються 
весь час незмінними. 

Введемо прямокутну декартову систему координат ОХ; х, Ха з початком у центрі Землі, осі Ох, 1 Ох, 


розташуємо довільно в екваторіальній її площині, а вісь Ох, спрямуємо за віссю її обертання. Для 
розпізнавання точок області у та її доповнення у перші будемо позначати грецькими, а другі - латин- 
ськими літерами, наприклад, 4 з 2) єу іха (зва у". Елемент об'єму замість 446,46, 
позначатимемо символом (С, а евклідову віддаль між точками хХ і 2 - у вигляді рез М Маси Землі 


М (2) С є у генерують у всьому просторі поле сили тяжіння, потенціал якого (при умові, що маси мають 
щільність ФМ (2) а с ме ), приймає вигляд 
ДС 
Ух) щ 7 пра ЗИє че 
У" В 5 5 


Ох), хєу 
1 
6 (1) 


, схеу"ї 
де 7 - гравітаційна стала, О(х) -0.5а) (х, Я х,) - потенціал ценробіжної сили, ш? - модуль вектора куто- 
вої швидкості Землі. Величина напруги поля або значення модуля градієнта потенціалу визначається так: 
3, ОМИ(х) дх, (х) 
шо ще Кк -/о У 
з(х)- УМ) - У, -(Ол)) 0) 
гі 23 Оп 
дх, (х ) я я - Й і 
де - оба со5(п, Хі 5 Ко -1,2,3 - напрямні косинуси одиничного вектора п (х) внутрішньої нормалі до 
Оп 
еквіпотенціальної поверхні «АУ(х.): У(у)- Сх, Сх з соп5і , яка проходить через точку Х, (зх) п (х) - 
скалярний добуток векторів 8! (х) іп (х). 
Зазначимо одну з важливіших характеристик модуля градієнта потенціалу, яка випливає з подання (2) і 
грунтується на наступному твердженні. 


Лема. Добуток двох функцій и(х) і »(х), хе класу С СД) буде гармонійною в області )) 
функцією ух) и(х) : у(х) тоді і лише тоді, коли кожний із співмножників буде гармонійним в Р, а їх 
градієнти Уи(х) і Ми(х) будуть ортогональними один до одного в області ПД). 

Доведення випливає негайно із рівності 
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УРу(х) б уз б у(х): У"и(х) В 2(Ми(х), Му(х)) ї и(х): У?м(х) 


Каі Ху 
Наслідок. Модудь градієнта потенціалу (сили тяжіння) не задовольняє рівнянню Лапласа в жодній 
точці області у 
ОМ (х) 
дх, 


у() ші со5(п, я ), К - 1,2,3 такими не є. У цьому легко переконатися безпосередньо. 


Насправді, хоч кожна із складових градієнта и(х) є гармонійною функцією, функції 


Якби на поверхні Землі (рівняння якої вважаємо заданим), крім значень модуля градієнта потенціалу 
«(х), хєду і внутрішньої нормалі т(х) до бу, вимірювався також напрямок градієнта п (о, то задачу 
відновлення потенціалу сили тяжіння можна було б звести до визначення потенціалу притягання 
У(х), хє у" у вигляді розв'язку зовнішньої задачі Неймана для рівняння Лапласа 








УРУ(х) - 0, хєу", п ш ФІ), хєду, У(Хх) -»00 (3) 
т 
де Ф(х)- о 2 сх )сов (п ої. со8(п, т) о Усі п, х, )сов(х., , т). 


Проте напрямок нормалі п (х) нам невідомий через виключну складність і високу вартість вимірю- 


вань, і від наведеного способу розв'язання задачі залишаються лише деякі міркування щодо побудови пев- 
ного ітераційного процесу. Перш ніж заглиблюватися в його деталі, охарактеризуємо коротко стан справ із 


граничними даними задачі відновлення потенціалу за значеннями модуля його градієнта 4(х) ш |- мУ(х) 5 


У рамках прийнятої моделі ці значення описуються, як неважко бачити, наступним виразом: 


коуч оМи ами 06 у 206 





Кеї дх, 
- гос) 
або 4 (х )- ах У, совіл,х, «чі з 
Кз! дх, | 
так як 90) о 
дх. 


Подамо потенціал притягання у вигляді суми УС)еу (х)-- Т(х) нормального (/ (х) і збурюючого 
Т(х) потенціалів. Нормальним будемо вважати потенціал, що генерується певними фіктивними масами, 
які дорівнюють масам в у , але на відміну від останніх розташовані в якомусь розумінні "правильно" (чи 
нормально?) в деякій іншій області у, З достатньо простою геометрією і з границею Оу, яка не дуже від- 
хиляється від земної поверхні Оу. У зв'язку з цим визначенням збурюючий потенціал буде характеризува- 
ти відхилення реального розподілу масв у від прийнятого нормального. Позначимо через У (х) одиничну 
(внутрішню) нормаль до поверхні С/ :/ (у)е С., а через У(х)е - Му (2) - модуль градієнта нормаль- 
ного потенціалу. Будемо вважати заданими напрямні косинуси со5(у, б) і со(х, ,т) внутрішніх норма- 

3 
лей У(х) та т(х) відповідно до поверхонь 0, 1 ду, аразом з ними і со5(у, т) - У сову, Хі )со8(х, 5 т) 
ке 

Тоді відновлення потенціалу притягання У(х), хе у" можна визначити на основі граничної задачі (3) 
у вигляді обчислення послідовних наближень УФО), Кк -- 0,1,2,...,00. Таким чином, за визначеними на 
попередньому і-му кроці наближеннями сок(п,, х, У К -1,2,3 напрямних косинусів нормалі л(х) визна- 
чимо на границі ду відповідно до формули (4) наступне і -- І-ше наближення сили притягання 

; чу 
АКРВИ ро РОЮ ВИН Ро АР 


па Ох, 





і значення функцій соб(п, у т) са У сов(п, Зк )соз(х, й т), хєду, 
геї 
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Фа (2 Чіа (х)сов(п, , т)- Ух) сову, т), хебу 


Далі розв'язуємо зовнішню задачу Неймана для рівняння Лапласа. 


Я ОТ, 
тан хе у ТК а хе т 


Її розв'язок будемо шукати у вигляді потенціалу простого шару мас, розподілених на поверхні ду з 





оо. (5) 


неперервною щільністю 6, (х), хєду, 


1 б. 
ла Є пе Мети хну" (6) 





Невідому щільність визначимо за граничною умовою, яка за відомою теоремою стрибка приводить 
до інтегрального рівняння шо се роду: 


РАКА Тк У, ел (2345, з Фа (2) хєду, (7) 


1 д ро бура бан 

2 дт, з 2л ро - 4 
ці 

бо 


У результаті розв'язку рівняння (7) обчислюємо з використанням співвідношення (6) наближення 
похідних потенціалу притягання 


ур(ху- ЛУ 20), СТ, а(х х) АХ п 1 х, -8і 5 


й дх, дх, дх, дх, 4я |с- «ї 





де ядро К (х, 2) з ,а соби, т), из-х- б визначається за формулою 





со. и, т) - Уа т, 35 








б. (2345., 2123 | (8) 


Із зовнішнього вигляду виразів (8) помічаємо, що похідні наближень збурювального потенціалу визна- 
чаються тільки у внутрішніх точні області у . У точках же її границі ду значення похідних у поданому 
вигляді визначити неможливо через наявність у підінтегральних функціях досить сильних, неінтегрованих 
особливостей. А між іншим, для продовження обчислень потрібно мати значення похідних збурювального 
потенціалу саме у точках границі ду. Щоб їх знайти, слід передбачити спеціальну регуляризацію інтегра- 
лів (8). Маючи на увазі, що така операція вже виконана і значення похідних (8) знайдені в точках хєду, 
продовжимо обчислення наступних наближень 

-1 
5 2272 
(пах А ОЮ), х є ду, чаа (2) вх) (пан) «у СО9СО д 
СО Кр ЧИ ПРЕСУ фо МО ЗХ АРК ЗР ГО 0 де СНОМ 
КІ дх 
Фо (2) -4Чіио (2) со5(ту » т) - у(х)сов(у, т), хе бу 


Визначивши ці наближення, знову перейдемо до розв'язку граничної задачі (5)-(7) зі зміною індексу 





і-- 1-го наближення збурюючого потенціалу Т.,, (х) ії наближення У ("") (х )- 7 (х )-Т, (5 ) потенціалу 


із 2 
притягання 1 т.д. 

Запропонована схема відновлення потенціалу притягання потребує точного обгрунтування. Необхідно 
навести переконливі міркування не тільки щодо розв'язку граничної задачі (6)-(7), якою замінена коректна 


. М а к - 
у математичному плані задача (5), але й довести збіжність самих наближень у! 102) до потенціалу притя- 
гання У(х), хе Б а також вказати метод обчислення похідних від наближень збурювального потенціалу 


у замкненій області у" 


Перш ніж перейти до дослідження окресленого кола питань, відмітимо, що теореми однозначного роз- 
в'язку задачі Неймана для рівняння Лапласа за допомогою потенціалу простого шару, встановлені на поча- 
тку століття, за модулем його градієнта зводяться головним чином до доведення збіжності обчислюваних 


наближень Ух), Кк - 0,1,2,...,со функції У(х), хєеу" 

Збіжність, наближень Ух) буде доведена, якщо пощастить виявити у свою чергу збіжність послі- 
довності т, (х)) наближень збурюючого потенціалу. Дійсно, якщо буде показано, що Т, (х) збігається до 
т(х) при К -»00, то звідси випливатиме не тільки збіжність УФ(х) до У(х), хє у", але і збіжність до 
своїх границь будь-яких інших наближень, що однозначно визначаються за Т, (х). Позначимо відношення 
.|Ут (« 
Мо(ої 


«У такому разі буде справедливою наступна теорема. 
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Теорема 1. Якщо я (х) - квадрат відношення модуля градієнта збурювального потенціалу до модуля 
градієнта нормального потенціалу - настільки мала величина порівняно з «(х) , що нею можна знехтувати, 
то послідовність розв'язків ПТ, (х)) граничних задач (5) збігається до збурювального потенціалу т(х) обла- 
сті у . 

Для доведення теореми досить показати, що збігається послідовність І б ху); тому що збіжність 


останньої, як це видно з подання (6), спричиняє збіжність послідовності т, (х)). А щоб установити збіж- 


І 


ність послідовності 5, (х)), відмітимо, що оператор Ад; тд, (х)-н Гк (ж, 2)5, (2)45 : рівняння (7) відпо- 
бу 








відно до теореми існування розв'язку мас обмежений обернений оператор А". Нехай |А" | «с«о0|, Ви- 
бі (х) - бі (У «2Ф,а (х) -Ф; (») 
урахуванням визначення функцій Ф, (х) може бути подана у вигляді 
5. 0)- 6.2 « 2с діл (2) - 4, (9. 
Покажімо, що першим доданком у правій частині нерівності можна знехтувати. У зв'язку з цим 
відмітимо, перш за все, рівність 
2 
з 20(х) ОТ, (х) от. (х) 2 
Я -1 у. ії 
4: (х)-У, ет зу" нд ччкнея м (МУ. 
ОХ дх, ду 
Якщо тепер вираз для градієнта функції Т,, (х) переписати з урахуванням очевидного співвідношення 


от, (х ) 


рої 71 (х)сов(п і у) УЗ УС), то прийдемо до наступної залежності 


ду 
«94 (ман ть | она 
аа У 

















ходячи з цього твердження, запишемо очевидну нерівність , Яка з 























діа (| сов(п з т) З со5(п ЯР: т) - 2сЇсоз(п і т) 








3 
У той же час із виразу со5(п, , У) 2, со5(п, ,, х, усо5(х,у)21- У" (х 
Ка! 


пленарне видно, що кут 


ду 





со(п, м) при виконанні умови теореми досить малий. Тому попереднє співвідношення з точністю до 


а х)- ТІ зу" «Дені і | (9) 


2 5 
Є (х) можна записати у вигляді 





ду 

2 

Тепер легко показати, що різниця сок(п,, т)-- сов(п, |,т) має порядок величини є (х). Насправді, 
оскільки за визначенням 1 з урахуванням наближеної рівності (9) маємо 


со5(п, ; т) і ор єна б щів о 5 т) - со5(у, т)- Ї -у? сеть | , (10) 


СО б, дх, ду 








5 2 Й й Р 
то, очевидно, згадана різниця за порядком не перевершує числа є (х). Далі, на основі подання (9) і з ура- 
хуванням граничної умови задачі (5) можна записати 
Чіа (х) 74; (2) - со5(и, та; (х)сов(п, | , т) "Чіа (х)сов(л, з , т)). 
З цього, якщо взяти до уваги умову теореми і рівність (10), одержуємо наступний ланцюжок нерівно- 
стей 


Ч: (х)- Чі (9 з 
Чіа (2)- Чі ОЇ з |соз(у, т) а М, (х)- Чо (9 : 
5. )- б, Со « 2сЇсов(у, т) й 4 (х)- чо (2, на основі якої пе- 


ічі 
реконуємося, що послідовність б Я (х)) швидко збігається у собі, якщо Їсох(у, т)| «1.У той же час, коли 





чу) 4, (9) є |оз(му ті) | 


«рози, т) 






































Так, остаточно маємо нерівність 


|сох(у, т) - 1, тобто, коли напрямок внутрішньої нормалі у(х) до еквіпотенціальної поверхні С/, збігає- 


ться (або йому протилежний) з напрямком нормалі т (х) до поверхні Землі бу, тоді, очевидно, наша зада- 
ча зводиться до добре вивченої зовнішньої задачі Неймана для рівняння Лапласа 1 може бути розв'язана за 


класичною процедурою. Виключаємо цей випадок із розгляду (як нецікавий і нереальний в даній ситуації) 
75 


Вісник Київського університету. Серія ,Геологія". Випуск 13. - 1995. 


1, вертаючись до послідовності і5 яі (ж), відмітимо, що із збіжності у собі елементарно встановлюється збі- 
жність до єдиної границі б (х) і тим завершуємо доведення. ш 

Проблема поширення формул (8), справедливих для внутрішніх точок області у , на граничні точки 
хєдбу знімається наступним твердженням. 


Теорема 2. Якщо щільність потенціалу простого шару (6) є неперервна на бу функція, то граничні 
значення частинних похідних першого порядку від потенціалу визначаються у вигляді 


ст(х) 1 рху ба 
Зх, і реа г 5(2)45, жо зб )соб(х, т, ) 


Доведення теореми елементарне і випливає негайно з подання потенціалу У вигляді 
І 45. б(х) 45 
т(х)а- --- Пд(2)- б 
нау ЛО вар Й 


ду | 
якщо диференціювання по х, тимчасово замінити диференціюванням по внутрішній нормалі т, до повер- 











хні ду у точці х, і прийняти до уваги відому тотожність 








180 1 ані 
Ї 45. з 0.5, хе у" 
дт, х-4| без 


Зауваження. Наведена формула визначення похідної погано пристосована для цієї мети через склад- 
ність земного рельєфу. Тому її і Й самого початку замінити еквівалентною формулою 


от (х Жуан 1 
а ЕР ее зб вк дя, або т). 
дх, ре-4| 2 

Індекси, якими у виразах (8) описано послідовність граничних задач, у обчисленні похідних опущені, 
як несуттєві у даному розгляді. 
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